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Abstract 

Background: Schistosomiasis transmission is typically focal. Understanding spatial variations of Schistosoma 
infections and their associated factors is important to help to invent site-specific intervention strategies. 

Methods: A five-year longitudinal study was carried out prospectively in 12 natural villages, Guichi district of Anhui 
province. A GIS-based spatial analysis was conducted to identify geographic distribution patterns of schistosomiasis 
infections at the household scale. 

Results: The results of the spatial autocorrelation analysis for 2005 showed that there were significant spatial 
clusters of human infections at the household level, and these results were in agreement with that of the spatial 
scan statistic. As prevalence of infections in humans decreased over the course of control, the spatial distribution of 
these infections became less heterogeneous. 

Conclusions: The findings imply that it may be necessary to re-assess risk factors of S. japonicum transmission over 
the course of control and to adjust accordingly control measures in the communities. 

Keywords: Schistosomiasis, Intervention, Geographical information system, Spatial-temporal distribution, 
Spatial autocorrelation analysis, Spatial scan statistic, Schistosoma japonicum 



Background 

Schistosomiasis, a snail-mediated helminthiasis, is one of 
the most prevalent parasitic diseases in the world [1]. The 
disease is endemic in approximately 70 countries with 
about 200 million people affected worldwide [2-4]. Docu- 
mented evidence indicates that Schistosoma japonicum 
has existed in China for over 2100 years [5]. Results of 
large-scale epidemiological surveys showed that schisto- 
somiasis was distributed in 12 provinces and an estimated 
11.6 million people were infected after the founding of the 
People's Republic of China [6]. Over the past 50 years sig- 
nificant achievements have been made on schistosomiasis 



* Correspondence: ybzhouiafudan.edu.cn; jiangqwiafudan.edu.cn 
'Department of Epidemiology, School of Public Health, Fudan University, 138 
Yixueyuan Road, Shanghai 200032, China 

2 Key Laboratory of Public Health Safety, Fudan University, Ministry of 
Education, Shanghai, China 

Full list of author information is available at the end of the article 



control through ongoing national control programs [7], 
To date, five out of 12 formerly S. japonicum-endemic 
provinces have reached the national criteria of transmis- 
sion interruption (i.e. no new infection), and the estimated 
number of human infections has been reduced from 
1,638,103 in 1989 (the first nationwide schistosomiasis 
sampling survey) to 726,112 in 2004 (the third nationwide 
schistosomiasis sampling survey), reflecting a 55.7% de- 
crease [8]. Despite these achievements, schistosomiasis is 
still of considerable economic and public health concern 
in China due to its potential for re-infection and recent 
reemergence in some previously controlled areas [5], In 
2004, two targets were established by the State Council of 
China - reduce the human prevalence of infection to less 
than 5% in all endemic counties by 2008, and then to < 1% 
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by 2015 [9]. In order to achieve these goals, an integrated 
control strategy, aimed at reducing the roles of bovines 
and humans as sources of infection for snails, has been 
conducted in a number of endemic regions of China [9]. 

Given the extensively ongoing national control pro- 
grams, an immediate question, which remains largely un- 
addressed, is how these control programs are changing 
local epidemiological characteristics of the disease and if 
the control programs are sufficient to suppress the trans- 
mission towards the set goals. To address the question, we 
need to first understand how local patterns of the disease 
change over the course of control; typically these patterns 
have both spatial and temporal dimensions. Advances in 
geographical information systems (GIS) and statistical 
methodologies for spatial analysis provide powerful tools 
to help characterize and improve our understanding of 
space-time distribution of diseases [10-16]. However, des- 
pite a large body of studies reporting spatial and/or tem- 
poral patterns of schistosomiasis infections, few studies 
have analyzed the space-time changes in the disease preva- 
lence at a finer scale (e.g., below county level), especially 
before and after an integrated national control strategy 
[7,10-12,17,18]. In this study, GIS-based spatial analyses 
involving spatial autocorrelation analysis and spatial scan 
statistic were employed to identify geographic distribution 
patterns of schistosomiasis at the household level. Further- 
more, changes in these patterns were evaluated over time 
in the Qiupu River region, where the integrated national 
control strategy was implemented. Knowledge on changes 
in the local epidemiological situation will aid in imple- 
menting and adjusting anti-schistosomiasis measures, and 
also increase the efficacy of these measures. 



Methods 

Study location 

The study took place in the Qiupu River region in 
Guichi district, Anhui province, China, which is, histor- 
ically, one of the most serious areas of S. japonicum in 
China (Figure 1). This study region has a population of 
approximately 32,000 people and 5.7 km of snail- 
infected area [19]. Most of the residents are farmers liv- 
ing along the Qiupu River, and their principal activity is 
the cultivation of rice. The residents frequently contact 
infested water due to their working and living circum- 
stances. Oncomelania hupensis, the intermediate snail 
host for S. japonicum, is distributed in the marshlands of 
the Qiupu River [19]. 

An integrated national control strategy was imple- 
mented in the Qiupu River watershed starting in 2006 
[9,19]. Details on the interventions were reported in our 
previous work [19] and the main measures included: 1) all 
bovines were removed from the study area in 2006; 2) all 
residents aged 5-65 years were checked annually for schis- 
tosomiasis infection (detected annually using indirect 
hemagglutination assay (IHA) and all IHA positives were 
examined using Kato-Katz thick smear test (three slides) 
to confirm infection, except for those who received direct 
fecal examination) and only the residents with positive egg 
counts were treated with a single oral dose of praziquantel 
(40 mg/kg body weight); 3) the sanitation facilities were 
improved by constructing lavatories and latrines and tap 
water was supplied; 4) a program of health education was 
implemented focusing on avoidance of snail-infected areas 
and the associated river water; 5) some snail habitats were 
managed by agricultural methods (e.g. plant wheat). 
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Cohort population 

Twelve natural villages (i.e. a hamlet, which is a small 
group of dwellings under an administrative village which 
typically consists of natural villages) were selected ran- 
domly from our study region in 2005 and followed until 
2009 (Figure 1). A cohort of population was selected 
from the 12 natural villages. The inclusion criteria for 
households were that: a) the household was from the se- 
lected natural village, b) the household resided in the 
natural village for more than 1 year, and c) the house- 
hold had residents aged 5-65 years. The inclusion criteria 
for residents were that an individual must be: a) a resident 
of the selected household, b) 5-65 years of age, and c) no 
plan to migrate out of the natural village during the study 
period. Within the cohort, each observed resident was 
georeferenced by household. Georeferencing was imple- 
mented using a hand-held global positioning system 
(GPS) (Mobile Mapper TM, Thales Navigation, Inc., 
USA). Each resident of the cohort had his/her stool 
(a fecal specimen each individual) examined annually dur- 
ing the autumn season using the microscopic, quantitative 
Kato-Katz thick smear (3 slides) technique [18], and the 
detail of this stool examination was reported by Zhou 
et al. 2011 [20]. The outcome measure was incidence for 
the defined cohort of population except for 2005 (a base- 
line survey, and the measure was infection prevalence). 
Residents with a positive fecal examination result were 
treated with a single dose of praziquantel (40 mg/kg). 

Written informed consent was obtained from all adult 
participants and from the parents or legal guardians of 
children aged < 18 years. Ethical approval for the study 
was endorsed by the Ethical Review Committee of School 
of Public Health, Fudan University. 

Data manipulation and statistical analysis 

Data were entered in Microsoft Excel version 2003 
(Microsoft Corp., Redmond, Wash.). A human S. 
japonicum infection was defined as the presence of at 
least one egg in three Kato-Katz smears. 

The cohort characteristics and infection prevalence/in- 
cidence were analyzed using the SPSS Statistical Package 
for Social Sciences (SPSS Inc., Chicago, IL, 2007). Hu- 
man prevalencewas analyzed using a generalized linear 
model (GLM) with a logit link and a binomial distribu- 
tion [21,22]. The base model included year (from 2005 
to 2009), sex, and age. 

Spatial autocorrelation analysis 

Both Global and Local spatial autocorrelation analyses 
were implemented in OpenGeoda 0.9.8.14 software 
(http://geodacenter.asu.edu/software). Spatial autocor- 
relation analysis was employed to detect significant differ- 
ences from a random spatial distribution of schistosomiasis. 
Moran's I spatial autocorrelation statistics were calculated 



and visualized in the form of Moran Scatter Plots. In order 
to take into account variance instability of infection preva- 
lences over time in both Global and Local spatial autocor- 
relation analyses, the empirical Bayes (EB) adjustment for 
the variance instability of prevalences was used [14]. The 
spatial weight files were created by using Distance-Band 
weight. A significance assessment was performed by means 
of a permutation test, and a reference distribution was gen- 
erated under the assumption that the prevalence was ran- 
domly distributed. The number of permutation tests was 
set to 9999 and the significant P-level was set as 0.05. 

Spatial cluster analysis 

Spatial analyses were performed in the SaTScan'" v8.2.1 
software (Kulldorff M. and Information Management 
Services, Inc.). Spatial scan statistics were used to detect 
and evaluate the clusters of schistosomiasis cases. These 
were done by gradually scanning a circle window across 
space and noting the number of observed and expected 
observations inside the window at each setting. The Ber- 
noulli model [15,16] for high prevalences was used, and 
the maximum spatial cluster size was set to 50% percent 
of the population at risk. Cases were defined as the 
population infected with 5. japonicum, and controls 
were the remainder of the population. Monte Carlo hy- 
pothesis testing [23] was used in this study and the 
number of replications was set to 9999 with significance 
at the 0.05 level. 

GIS mapping 

A digitized polygon map was obtained for the Qiupu River 
watershed at a scale of 1:50,000. Clusters were mapped 
using ArcGIS 9.2 software (Environmental Systems 
Research Institute, Inc., Redlands, CA, USA) in order to 
identify their physical location. 

Results 

Cohort characteristics 

The population characteristics of the observed cohort 
are shown in Table 1. The mean age of the population 
increased by year, except for 2009 (F = 6.15, P =0.000). 

Infection prevalence and global spatial autocorrelation 

The infection prevalence of the cohort population de- 
creased from 9.3% in 2005 to 3.7% in 2009, reflecting a 
60.2% reduction. By mean of GLM test, this reduction 
was clear after adjusting for age and gender (Wald x = 
37.81, P = 0.000). The positive individual was found in 
every natural village in both 2005 and 2008. No infec- 
tions were found in 1 (in both 2006 and 2009) or 2 (in 
2007) natural villages in the other years. There was only 
one natural village where no positive individuals were 
found during two surveys. The Moran's I value for 
households decreased from 0.083 in 2005 to -0.022 in 
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Table 1 Characteristics of the cohort population examined from 2005 to 2009 



Year Number of Mean of households per Number of Mean of population per Mean of residents per Mean age 

households natural village population examined natural village household (year) 



2005 


275 


22.9 
(16-35) 


654 


55 
(33-85) 


2.38 


374 
(6.7-65.0) 
[36.2-38.6] 


2006 


277 


23.1 
(16-35) 


652 


54 
(33-84) 


2.35 


38.2 
(7.7-645) 
[37.1-39.4] 


2007 


275 


22.9 
(16-35) 


649 


54 
(33-83) 


2.36 


40.2 
(9.7-65.0) 
[39.0-41.4] 


2008 


261 


21.8 
(12-32) 


575 


'18 
(25-85) 


2.20 


41.2 
(7.3-65.0) 
[39.9-42.4] 


2009 


268 


22.3 
(13-35) 


684 


5/ 
(27-85) 


2.55 


38.2 
(5.9-65.0) 
[36.9-39.4] 



The range min-max values in parentheses (), 95% CI in square bracket []. 



2009. There was a significant global spatial autocorrelation 
only in 2005 (P<0.05, 9999 permutations). No significant 
global spatial autocorrelation was observed in the other 
years (P > 0.05, 9999 permutations). These results indi- 
cated that the spatial distribution of S. japonicum infec- 
tions might change from correlation to dispersion. The 
trend of the Moran's I value at the household level is very 
similar to that of the schistosoma infection prevalence 
(Figure 2). 



Household (Moran's I) 



Prevalence/1 ncidence 



0.10 



0.05 



0.10 



0.00 




-0.05 



Figure 2 Global spatial autocorrelation analyses for annual 
prevalence/incidence of infection S. japonicum from 2005 to 
2009. Prevalence/Incidence represents that of the cohort 
population. The first point (2005) in Prevalence/Incidence represents 
infection prevalence of the cohort, and the remaining points 
represent infection incidence of the cohort. The points in Household 
(Moran's I) represent the Moran's I values at household level. 



Local spatial autocorrelation 

Based on the values of the global Moran's I above (i.e. 
the global Moran's I value (0.0376) of 2006 is very close 
to that (0.0426) of 2007 and the global Moran's I value 
(-0.0246) of 2008 is very close to that (-0.0219) of 2009) 
(Figure 2), the data from 2005, 2007, and 2009 were 
selected to analyze the local spatial autocorrelation of 
5. japonicum transmission. Local Indicators of Spatial As- 
sociation (LISA) Moran's I spatial autocorrelation analyses 
for S. japonicum infection prevalence are shown in 
Figure 3, in which the horizontal axis is the standardized 
value for S. japonicum infection prevalence, and the ver- 
tical axis is the mean standardized neighbor value 
(W RATE). These Moran scatter plots provide a visual 
representation of spatial autocorrelation. These plots were 
divided into four quadrants to allow a qualitative assess- 
ment of the association by type: high-high (upper right) 
and low-low (lower left) as positive local spatial autocor- 
relation and spatial clusters, and high-low (lower right) 
and low-high (upper left) as negative correlation and 
spatial outliers. 

Most of the households observed were distributed in 
the low-high quadrant for the three years examined, and 
the number of the households in high-high quadrant 
decreased each year, indicating the spatial distribution of 
S. japonicum infections became less clustered with time 
(Figure 3). The proportion of significant clusters or out- 
liers in the observed households decreased each year 
(P<0.05), except for the high-low outliers (j 2 =1.57, 
P>0.05). No significant high-high clusters were ob- 
served in 2009 (Table 2). 
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Figure 3 LISA Moran scatter plot for S. japonicum infection prevalence at the household level: (a) 2005, (b) 2007, (c) 2009. 



Space cluster analysis 

Only one significant spatial cluster of infection was ob- 
served in 2005 (P < 0.01) (Figure 4), no significant spatial 
clusters were found in the other years (P>0.05). The 
GPS coordinates of the center of the significant spatial 
cluster was 30.57686 N and 117.37817 E, and the radius 
of the cluster was 1.44 Km. This cluster included 36 
cases and the infection prevalence of the population in- 
side this cluster was 18.5%, while outside the cluster, 
there were only 25 cases with 5.4% infection prevalence. 
The cluster was near the Qiupu River (Figure 4). 

Discussion 

Transmission of S. japonicum is influenced by many 
factors. Some key factors are climatic suitability (e.g. 
temperature and precipitation); spatial and temporal dis- 
tribution of the intermediate host, O. hupensis; human 
activities (e.g. occupational, domestic, and recreational 
water-contact patterns); environmental contamination with 
human and animal excreta, in particular, cattle stool; and 
lack of a clean water supply, sanitation, and proper hy- 
giene [12,24-28]. As a result, S. japonicum infection 

Table 2 Proportion of significant clusters or outliers in 
the observed households at a = 0.05 level, 9999 
permutations 



Scatter plot 




Year (%) 




X 2 


P 


quadrant 


2005 


2007 


2009 




value 


High-high 


2.55 (7/275) 


0.73 (2/275) 


0.00 
(0/268) 


8.81 


0.013 


High-low 


1 .09 (3/275) 


0.36 (1/275) 


0.37 
(1/268) 


1.57 


0.456 


Low-low 


19.27 (53/275) 


10.55 (29/275) 


1.49 
(4/268) 


45.61 


0.000 


Low-high 


8.36 (23/275) 


7.27 (20/275) 


3.36 
(9/268) 


6.30 


0.043 



exhibits marked spatial heterogeneity from the community 
scale even at the scale of a single administrative village) up 
to the regional scale, [7,17,28-31]. Our results from the 
spatial autocorrelation analysis showed that there was sig- 
nificant spatial cluster at the household level in 2005, 
when the infection prevalence was relatively high (i.e. > 
9.3%). Similar results were obtained from analyses of the 
spatial scan statistic that was used to identify areas and 
population with high prevalence (a significant spatial clus- 
ter of infection was found at the household level in 2005). 
These findings were similar to the results reported by Peng 
et al. (2010) [17], where the human infection prevalence 
was 6.5%. Additionally, the results support the general view 
of the focal nature of S. japonicum (i.e. S. japonicum in- 
fection shows distinct spatial distribution (focality and 
aggregation) when the human infection prevalence is rela- 
tively high). 

However, the four years' results of both the global 
spatial autocorrelation analysis and the spatial scan 



The number of significant clusters (or outliers) and observed households are 
presented in parentheses. 




A 

Water body 



Figure 4 Spatial cluster location of S. japonicum transmission at 
the household level in 2005 year. Black dots represent the 
households. Rate Ratio is presented in the parenthesis. 
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statistical analysis all showed that no significant spatial 
cluster of infection was found at the household level 
after 2005, the results of the local spatial autocorrelation 
analysis also showed that the spatial distribution of 
S. japonicum infections became less clustered with 
time. These suggested that the spatial distribution pattern 
of S. japonicum infections might be changing over the 
course of control, from a heterogeneous pattern to a more 
homogenous pattern with a decrease in the infection 
prevalence at the household scale in our study area. This 
was different from the general phenomenon of infectious 
diseases (i.e. spatial heterogeneity increases with declining 
rates of transmission). This difference might be attributed 
to or related to 1) the complex life cycle of S. japonicum 
involving snail (Oncomelania hupensis) intermediate hosts 
and many vertebrate definitive hosts [24,32] (for example, 
the geographic distribution of schistosomiasis japonica 
largely depends on the availability of the susceptible 
schistosome- transmitting snail O.hupensis [24], and the 
spatial distribution of snails infected with S. japonicum be- 
comes less heterogeneous with a decrease in the infection 
prevalence [33]); 2) the effective integrated national inter- 
vention being employed during the study period [9]; and 
3) changes in some important risk factors influencing the 
transmission of S. japonicum in our study area. Although 
some control measures of snails were carried out in our 
study region, the spatial distribution of snails does not ob- 
viously change, in other words, the area with snails does 
not change over the course of control [19]. However, the 
integrated national control strategy resulted in changes in 
some important risk factors influencing the transmission 
of S. japonicum in the study area [19]. For example, envir- 
onmental contamination with cattle excreta was com- 
pletely wiped out after all bovines were removed from the 
study region. Previously, buffaloes were responsible for ap- 
proximately 75% of human transmission [34]. Hygiene 
was improved, to a large extent, by constructing lavatories 
and latrines and by supplying tap water. Furthermore, new 
human infections of schistosomiaisis each year were ef- 
fectively cleared by treatment with praziquantel. Also, 
health education programs might have changed human 
behavior. However, it is unclear whether there would be 
changes in other risk factors with the shift from clustered 
to random spatial distribution. For example, whether 
some minor and neglected risk factors (e.g. environmen- 
tal contamination with wild (e.g. mouse and rabbit) and 
other domesticated (e.g. pig and dog) animals excreta and 
population movement) would become important risk 
factors of S. japonicum transmission when other key 
risk factors noted above were controlled. This is neces- 
sary to further survey and to reassess the risk factors of 
S. japonicum transmission. The findings of random 
spatial distribution patterns of S. japonicum transmission 
showed that it might be necessary to adjust current 



control measures in the community due to the changes in 
risk factors of S. japonicum transmission. For example, if 
some neglected risk factors (e.g. wild mouse, population 
movement) become important risk factorsfor 5. japonicum 
transmission in a community after bovines are removed as 
important sources of infection, an approach for wild rat 
control or systematic administration of antischistosomal 
drugs (e.g. praziquantel) to migrantsmay need to be con- 
ducted in the community. This can help to tailor effective, 
locally adapted control measures [27] . 

Utzinger et al. (2003) [27] reported that the spatial dis- 
tribution of S. mansoni infection intensity levels was ran- 
dom or homogenous in a single village, but the infection 
prevalence of S. mansoni was very high (80.4%) when 
compared with the infection prevalences of S. japonicum 
in our study. In our study, we saw a random spatial pat- 
tern of S. japonicum infection only when the human in- 
fection prevalence was relatively low (< 4%); when the 
infection prevalence was relatively high (e.g. > 9%), the 
spatial pattern was heterogeneous. These findings indi- 
cate that spatial patterns of Schistosoma transmission 
might be relevant to infection prevalence at a finer scale, 
but it is unclear whether there is a critical prevalence 
when the distribution shifts from heterogenous to ho- 
mogenous distribution. Additionally, these results sug- 
gest that differences in spatial distribution patterns of 
schistosomiasis might be related to differences in the 
spatial distribution of risk factors of Schistosoma trans- 
mission, this needs further study. 

The cohort population was a stable population except 
in 2008 when 8.8% of the participants were lost to 
follow-up. This was mainly due to population movement 
(e.g. mobility for paid work or school). There was no ob- 
vious difference (e.g. age, sex, occupation) between the 
observed population and the population lost to follow- 
up (data not shown). A recognized limitation of this 
study was related to the Kato-Katz thick smear tech- 
nique used. This method has become relatively insensi- 
tive due to widespread chemotherapy that results in 
generally low worm burdens [35]; however, our previous 
study [36] showed that the sensitivity of a single stool 
examination with the three-slide Kato-Katz method did 
not change significantly when the prevalence of a single 
stool examination with the three-slide Kato-Katz method 
is between 13.0% (close to the baseline (9.3%) in 2005) 
and 3.9% (close to the prevalence of 3.7% in 2009). Sig- 
nificant spatial clusters of infection could be found by 
both the spatial autocorrelation analysis and the spatial 
scan statistical analysis in 2005; however, significant 
spatial clusters of infection could not be found after 
2005. So, the relative insensitivity of Kato-Katz method 
might not be a major factor for the apparent lack of 
clustering. Other weaknesses of the study include a rela- 
tively small dataset and some risk factor data related to 
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households or individuals were not collected (e.g. human 
activities and hygiene). This limits the ability to further 
analyze why the spatial distribution of 5. japonicum infec- 
tions changed. Given the results of the present study, it is 
cautiously suggested that, at the household level, the 
spatial distribution of S. japonicum infections become less 
heterogeneous as the prevalence of infection decreases. 

Conclusions 

The spatial distribution of 5. japonicum infections might 
become less heterogeneous as the prevalence of infection 
decreases at a household scale in our study field. The risk 
factors of S. japonicum transmission might change in our 
study field. These imply that it may be necessary to re- 
assess risk factors of 5. japonicum transmission over the 
course of control and to adjust control measures accord- 
ingly in the communities. 
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